CASE FILE 
CQpy 



NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 



REPORT No. 636 



APPROXIMATE STRESS ANALYSIS OF MULTI- 
STRINGER BEAMS WITH SHEAR DEFORMATION 

OF THE FLANGES 



By PAUL KUHN 



REPRODUCED BY 

NATIONAL TECHNICAL 
INFORMATION SERVICE 

U.S. DEPARTMENT OF COMMERCt 
SPRINGFIELD, VA. 22161 



1938 




AERONAUTIC SYMBOLS 
1. FUNDAMENTAL AND DERIVED UNITS 



Length 
Time__ 
Force.. 

Power _ 
Speed. 



Symbol 



I 
t 

F 



P 
V 



Metric 



Unit 



meter 

second 

weight of 1 kilogram _ 



horsepower (metric) . 
fkilometers per hour. 
\meters per second _ _ 



Abbrevia- 
tion 



m 

s 

kg 



k.p.h. 
m.p.s. 



English 



Unit 



foot (or mile) 

second (or hour). _ 
weight of 1 pound 

horsepower 

miles per hour 

feet per second — 



Abbrevia- 
tion 



ft. (or mi.) 
sec. (or nr.) 
lb. 



hp. 

m.p.h. 
f.p.s. 



2. GENERAL SYMBOLS 



m, 

i 



s, 
o, 

b, 
c, 
b 2 
S' 
V, 

2. 

L, 

D, 

D 0 , 

Do 

D„ 

C, 

B, 



Weigh t—mg 

Standard acceleration of gravity = 9.80665 
m/s 2 or 32.1740 ft. /sec. 2 
W 

Mas3= — 
g 

Moment of inertia=ra£ 2 . (Indicate axis of 

radius of gyration k by proper subscript.) 
Coefficient of yiscosity 

3. AERODYNAMIC SYMBOLS 



v, Kinematic viscosity 
p, Density (mass per unit volume) 
Standard density of dry air, 0.12497 kg-nv*-s* at 

15° C. and 760 nun; or 0.002378 lb.-ft." 4 sec. 2 
Specific weight of "standard" air, 1.2255 kg/m 3 or 
0.07651 lb./cu. ft. 



Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio 

True air speed 



Dynamic pressure —i^pV 2 

Lift, absolute coefficient ^L^~g 



Do 



Drag, absolute coefficient @d = ~£ 
Profile drag, absolute coefficient C^o^^g 
Induced drag, absolute coefficient Cz>i— ~g 
Parasite drag, absolute coefficient C Dp =^ 



C 



Cross-wind force, absolute coefficient Cc=jjjg 
Resultant force 



ft* 

Q, 



VI 



7, 



Angle of setting of wings (relative to thrust 
line) 

Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 

Reynolds Number, where I is a linear dimension 
(e.g., for a model airfoil 3 in. chord, 100 
m.p.h. normal pressure at 15° C, the cor- 
responding number is 234,000; or for a model 
of 10 cm chord, 40 m.p.s., the corresponding 
number is 274,000) 

Center-of-pressure coefficient (ratio of distance 
of c.p. from leading edge to chord length) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) 
Flight-path angle 



REPORT No. 636 



APPROXIMATE STRESS ANALYSIS OF MULTI- 
STRINGER BEAMS WITH SHEAR DEFORMATION 

OF THE FLANGES 

By PAUL KUHN 
Langley Memorial Aeronautical Laboratory 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



HEADQUARTERS, NAVY BUILDING, WASHINGTON, D. C. 
LABORATORIES, LANGLEY FIELD, VA. 

Created by act of Congress approved March 3, 1915, for the supervision and direction of the scientific 
study of the problems of flight (U. S. Code, Title 50, Sec. 151). Its membership was increased to 15 by 
act approved March 2, 1929. The members are appointed by the President, and serve as such without 
compensation. 

Joseph k\ Ames, Ph. D., Chairman, 

Baltimore, Md. 
David W. Taylor, D. Eng., Vice Chairman, 

Washington, D, C. 
Willis Ray Gbegg, Sc. D., Chairman, Executive Committee, 

Chief, United States Weather Bureau. 
William P. MacCbaoken, J. D., Vice Chairman, Wxecutiwi 
Committee, 

Washington, D. C. 
Cb m;li s (;. Abbot, Sc. i >. 

Secretary, Smithsonian Institution. 
Lyman J. BbIQGS, Ph. D., 

Director, National Bureau ol Standards. 
Arthur B. Cook. Rear Admiral, United States Navy, 

Chief, Bureau of Aeronautics, Navy Department. 
iiakky F. Guggenheim, m. A., 

Porl Washington, Long [sland, N. V. 



George W. Lewis, Director of Aeronautical Research 
John F. Victory, Secretary 
J. E. Reid, Kngineer-in-Charge, Langley Memorial Aeronautical Laboratory, Langley Field, Va. 
John J. Ide, Technical Assistant in Europe, Paris, France 



Sydney M. Kraus, Captain, United States Navy, 

Bureau of Aeronautics, Navy Department. 
Charles A. Lindbergh, LL. D., 

New York City. 
I n \ is Mulligan, J. S. 1 >.. 

Director of Air Commerce, Department of Commerce. 
Augustine W. Robins, Brigadier General, United States Army. 
Chief Materiel Division, Air Corps, Wright Field. 

I >;iYton, ( )l)io. 
EDWARD P. WabNEB, M. S.. 

Greenwich, Conn. 
Oscab Westoveb, Major General, United States Army, 
chief of Air Corps. War Department. 

( >BVTLLE WBIGHT, Sc. D., 

Dayton, Ohio. 



TECHNICAL COMMITTERS 



aerodynamics 

POWER PLANTS FOR AIRCRAFT 
AIRCRAFT MATERIALS 



AIRCRAFT STRUCTURES 
AIRCRAFT ACCIDENTS 
INVENTIONS AND DESIGNS 



Coordination of Research Needs of Military and Civil Aviatioyi 
Preparation of Research Programs 
Allocation of Problems 
Prevention of Duplication 
Consideration of Inventions 



LANGLEY MEMORIAL AERONAUTICAL LABORATORY 

LANGLEY FIELD, VA. 

Unified conduct, for all agencies, of 
scientific research on the fundamental 
problems of flight. 



OFFICE OF AERONAUTICAL INTELLIGENCE 

WASHINGTON, D. C. 

Collection, classification, compilation, 
and dissemination of scientific and tech- 
nical information on aeronautics. 



REPORT No. 636 



APPROXIMATE STRESS ANALYSIS OF MULTISTRINGER BEAMS WITH SHEAR 

DEFORMATION OF THE FLANGES 



By Paul Kuhn 



SUMMARY 

The problem of skin-stringer combinations used as 
a.rlalhj loaded jut m Is or as con rs for box bra //is is con- 
sidered from the point of rieto of the practical stress 
analyst. By a simple substitution the problem is reduced 
to the problem of the single-stringer structure, which 
has been treated in N. A. C. A. He port Ao. (JOS. The 
method of making this substitution is essentially empiri- 
cal; in order to justify it, comparisons are shown between 
calculations and strain-gage tests of three beams tested 
by the author and of one compression panel and three 
beams tested and reported elsewhere. 

INTRODUCTION 

A combination of a plate and stringers is frequently 
used as a structural element. Figure 1 (a) shows such 
a combination used as a tension member; figure 1 (b) 
shows one used as the tension side of a beam. The 
stress distribution in structures of this type is materi- 
ally influenced by the shear deformation of the plate. 
In aeronautical structures, where the plate often con- 
sists of a thin sheet that may be allowed to buckle 
into a diagonal-tension field, it becomes necessary to 
consider the effect of this shear deformation more 
carefully than is customary in other types of structure. 

Reference 1 discusses in detail (lie fundamental prin- 
ciples and the simplifying assumptions that permit a 
mathematical approach to the solution of the problem. 



(a) 

Figure 1.- 




Skin-st ringer combinations as structural elements. 



It is shown that numerical solutions can be obtained if 
there is only a single central stringer (fig. 2). A 
thorough familiarity with the method of analyzing 
single-stringer structures as given therein is presup- 
posed. For multistringer structures the mathematics 
becomes so complex that there is very slight possibility 
of obtaining sufficiently general solutions on the basis 



of the assumptions thai were used for the single- 
si ringer si rucl ures. 

Methods combining a desirable degree of accuracy 
with a reasonable degree of generality will, in all 



^^^^^^^ 




(b) 

Figure 2.— Single-stringer structures. 



probability, be methods of successive approximation. 
Attempts to develop such a method have thus far 
failed because the convergence is prohibitively slow . 
When such a method is found, it is not likely to be very 
rapid. Approximate methods developed in the interim, 
such as the one to be presented in this paper, will 
therefore retain their value by furnishing a very useful 
first approximation. 

The method presented herein w r as devised to answer 
the urgent need for estimating the effects of shear 
deformation. It aims chiefly at rapidity and ease of 
application, which are achieved at the expense of intro- 
ducing some empiricism. The experimental evidence 
presented is believed to be sufficient to prove that the 
method depicts reasonably well the influence of the 
shear deformation on the stringer stresses. 

METHOD OF ANALYSIS 

It is customary to designate tensile stresses and forces 
as positive. Figures, derivations, and formulas pre- 
sented herein deal, in general, with tension members. 
The only differences between tension members and 
compression members are quantitative differences in 
the effective widths and in the effective shearing stiff- 
nesses of the sheet. In the case of beams, the side not 
under consideration at the moment, i. e., the compres- 
sion side in most of the discussion of this paper, is 
assumed to be concentrated at the shear web (figs. 
1 (b) and 2 (b)) in such a location that the effective 
depi \\ is not changed. 
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The investigation of reference 1 was restricted to 
symmetrical structures as indicated in figures 1 and 2. 
The same restriction will be made in the present paper, 
and formulas and numerical data must be understood 
to apply to the half structure unless otherwise specified. 

In order to unify the terminology, the designations 
and symbols used in reference 1 for beams are extended 
in the present paper to axially loaded panels. (Sec 
appendix A for a list of symbols.) The directly loaded 
stringer of an axially loaded panel will therefore be 
referred to as the "flange" (subscript F) and the other 
stringers attached to the sheet as "longitudinals" or 
"stringers" (subscript L). This procedure is justified 
because the axially Loaded panel may be considered as 
the cover of a box beam in pure bending under the 
assumptions made. 

It is assumed in all cases that the longitudinals are 
distributed uniformly along the chord. It is furthermore 
assumed that camber is moderate, not exceeding the 
amounts found, for instance, in wing beams. Finally, it 
is assumed that the effective shear stiffness and the sheet 
thickness are constant along the chord. 

GENERAL PRINCIPLES OF METHOD OF ANALYSIS 

The mathematics of the mult isl ringer beam with 
\ ariable cross section is too complex to admit of ready 
solution. Broadly speaking, two methods of procedure 
may be used in such a case. One method would be to 
use approximate methods of solving the equations; the 
other method would be to ideali/e and simplify the 
physical concept of the structure until the mathematical 
relations become manageable. The second method is 
used in this paper. 

The results obtained in reference 1 show that the 
highest stresses occur at the flange and that they 
decrease from the flange toward the center line of the 
structure. The stress in the flange and the closely 
related stress in the longitudinal adjacent to the flange 
are therefore of paramount interest to the analyst. 

In beams with cambered cover, which were not 
treated in reference 1, the highest stress in the longi- 
tudinals may occur adjacent to the flange or it may 
occur at the center line of the beam. When it occurs 
at the center line, the stress there also becomes a 
matter of concern to the analyst. 

It is quite obvious that, in general, the most impor- 
tant physical actions will take place around the flanges, 
partly because the loads are applied there and partly 
because the stresses reach a maximum there as long as 
there is no violation of the basic requirement that the 
camber be very moderate. Consequently, any sim- 
plification that may be made should affect as little as 
possible the picture of the physical relations in the 
immediate vicinity of the flanges. 

In conformance with this requirement, the simplifica- 
tion necessary for obtaining a solution was achieved by 
using a "substitute structure" obtained by leaving the 



flange {and shear web) intact but replacing the longitu- 
dinals that are actually uniformly distributed over the 
width of the sheet by a single longitudinal equivalent to 
them as far as action on the flange is concerned. This 
substitution reduces the problem of the multistringer 
structure to that of the single-stringer structure, which 
can be analyzed as shown in reference 1. The method 
of substituting (temporarily ) a simplified structure for 
the actual one corresponds in part to the method of 
using "phantom members" in trusses. 

The substitute structure is used only to calculate the 
stresses in the part that it has in common with the 
actual structure, namely, the flange and the skin ad- 
jacent to the flange. After this object has been 
attained, the substitute structure is discarded. The 
stresses in the actual distributed longitudinals are then 
obtained by using the method described in reference 1 
for distributing "corrected forces. " 

It is clear that, in any given case, at least one equiva- 
lent single longitudinal exists. Whether or not there is 
a general method for finding this equivalent longitu- 
dinal, however, is a question that could be answered 
theoretically only if all the exact mathematical solu- 
tions were known. They are not known, and the 
method of finding the equivalent longitudinal is there- 
fore essentially empirical and must be justified by tests. 
This requirement is not such a serious drawback as it 
may seem to be, because the basic simplifications used 
are such that experimental verification is required many 
event. 

The method of finding the equivalent single longitu- 
dinal isasfollows: Remove from the sheet each individual 
stringer of cross-sectional area A at a distance y from 
the center line and attach, at the center line of the 
structure, a substitute stringer with a cross-sectional 
area 

A S =A^ (1) 

w here a y is the stress in the actual stringer and v C l the 
stress in the actual center-line stringer. The ratio 
GuI°cl may be considered as the "effectiveness'' of the 
stringer at >/ relative to the stringer at the center lint 4 
y=0; the use of this factor in expression (1) tends to 
counteract the loss of effectiveness caused by moving 
the stringer from its original location to the center line. 
The sum of the individual substitute stringers attached 
at the center line constitutes the single equivalent 
longitudinal. 

As the stresses a, and a CL are unknown at the outset, 
for a first approximation, the ratio g v \gcl is obtained 
from equation (17) of the constant-stress solution given 
in reference 1. With the stresses thus computed, a 
second approximation might be made. In all cases in- 
vestigated thus far, it was found that the second approx- 
imation agreed with the first one within the limits of 
experimental accuracy. The use of the second approxi- 
mation is therefore considered unnecessary. (It must 
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be borne in mind that the method of finding the equiva- 
lent longitudinal is essentially empirical. Conse- 
quently, there is no valid reason to believe that the 
second approximation must be better than the first one.) 

ANALYSIS OF AXIALLY LOADED PAN 

As an example of the analysis of an axially loaded 
panel, the una lysis of the compression panel with seven 
si if! mors, described in reference will he discussed in 
detail. The pertinent data on this panel are given in 
figures 3 (a) and 3 (b). 

Estimate of effective areas and of effective shear 
stiffness. — The test results are given in reference 2 for 
•2/' = 2,000, 4,000, and 6,000 pounds. The analysis will 
be made for 2P=4,000 or P = 2,000 pounds. It will 
become apparent that the conditions at this load are 
the same as for very small loads, so that the analysis 
will be valid for any load between 0 and 4,000 pounds. 

The mean stress in the panel (reference 2) is 



2,000 
" 0.70 



= 2,860 lb./sq. in. 



This stress is fairly close to the compressive buckling 
stress of the sheet; the effective width of the sheet will 
therefore be taken as equal to the actual width. The 
effective stringer area for the llange is therefore 

^==0.180+2X0.024 = 0.228 sq. in. 

and for the sum of the other stringers 

Al=2.5 X 0.088 + 1 0 X 0.024 =0.460 sq. in. 

The force at the bottom of the edge stringer is approxi- 
mately 

7<V=2,860X0.228=652 lb. 

leaving 1,348 pounds to be transmitted by shear in the 
sheet to the other stringers. The average shear stress 
in the sheel next to the edge stringer is therefore 

T= 48xS(J24 = 1 ' 17() lb,/sq - m ' 

The critical buckling stress for 0.024-inch dura! sheet, 
t inches wide and assumed simply supported, is, accord- 
ing to Timoshenko, 

r C ri*=l,730 lb./sq. in. 

This value is so far above the actual stress that there is 
no possibility of a diagonal-tension field forming and 
reducing the shear stiffness, so that QJE—QA0 may be 
taken. 

Determination of substitute structure. — Figure 3 (c) 
shows the cross section of the idealized structure 



assumed for the analysis. The stringer areas given on 
this figure are effective areas that include the effective 
width of the sheet; the sheet is now assumed to carry 
only shear. 

If there are at least two intermediate stringers be- 
tween the center stringer and the llange, the calculation 
of the substitute stringer may be simplified by using a 
formula derived on the assumption that there are in- 
finitely many intermediate stringers; that is, on the 



P= 2,000 I Symmetrical about <£ 

r 



t = 0.024 



W////« 



A' 0.180 



A- 0.088 



h£ — 4^ — 4 — H 1 



• 4 + 4 — >| — 4 \ 

(b) t 



L =48 





/ *0.02* 


t 


-0.226 A ^ 0460 


• » 




t r - 




h 


< — 4 >jc — 4 > 


c — 4~ 





(c) 



<-b = /2-> 
(a) 



A F = 0.228 



-t- 0.024 



A LS = 0.500^ 



-f2 - 
(d) 



Figure 3.— Compression panel used for sample analysis. 

assumption that the area A L of the intermediate 
stringers is distributed uniformly along the width b of 
t he sheet. The derivation of this formula is ns follows: 
According to the constant-stress solution (reference 1, 
equation (17)) 



— = cosh K 3 y 



(2) 



where 



and, in the case of a constant cross section, 
KJj 



2A L bE 



tuo, 



(3) 



The area of an individual stringer is now 

and the area of the substitute stringer that replaces 

it al tlx* center line is, according to equation (1), 
dAzs^-ydy cosh K s y 



4 



REPORT NO. 636 NATIONAL ADVISOR!? COMMITTEE FOB AKRONAUTK's 



The total area of the substitute stringer located at the 
center line is therefore 4 

A LS = ~yJ^ cosh K 3 y <fy=A L K ^ (4) 

In the case under consideration 

m 

so that 



I 2X0-460X12 =0>706 
\ 0.024 X48 2 X 0.40 



^=0.460X^^=^0.500 sq. in. 

Figure 3 (d) shows the cross section of the substitute 
structure. 

Analysis of substitute structure. — The substitute 
st ructure of figure 3 (d) can be analyzed by applying the 
Formulas given in appendix B. By formula (A-l) 

_0.40X0.024/ 1 1 \ 

12 \0.228"*" 0.500/ 

£l=0.0715 



rtAx 



x+Ax 



Figure 4.— Free-body diagram for calculating shear stress. 

For any station along the span, the stresses and forces 
can now be calculated. For example, ai the bottom of 
the panel O 0), by formula (A 3) 

P A , A™ cosh K 



I / 1 ; > COStl h r v 
**~~A w -\ ; A, bosh KL) 



2,000 / 0.50 
.500V + 0.22<< 



>oo 

28 • 



[.00 \ 
5.53 ' 



0.228H 0. 
=3,134 lb./sq. in. 

With the computation of o> the substitute structure lias 
served its purpose and is discarded. It is important 
not to confuse it with the actual structure in any of the 
following computations. 

Calculation of stresses in longitudinals. — The total 
force F L in the actual longitudinals is 

F l = P-F f =P-g f A f 

or, at x=0, 

Fi=2,000— 3,134X0428=1,286 Lb. 

This force is to be distributed over the longitudinals 
by the method given in reference 1, equations (21) and 



(22). In order to apply this method, compute the 
average stress 



=^=S=2,8001b./sq. in. 



A L 0.460 



and the ratio 



Q-L ag ^ 2,800 
<r F 3,134 = 



0.894 



With this ratio as abscissa, read from figure IS of ref- 
erence 1 (redrawn to a larger scale) 



76=0.605 



and calculate 



a> _3,134_ 



=2,040 lb./sq. in. 



* CL cosh Yb 1.188 = 



For the other two stringers, which are located 
at y=Ysb and y=%b, the stress will be, for stringer C, 

a=a CL cosh Fy=2,640Xcosh 0.202 = 2,694 lb./sq. in. 

and, for st ringer B, 

cr = 2,640Xcosh 0.404 = 2,860 lb./sq. in. 

The shear stress r at any point in the sheet is obtained 
most conveniently by considering the equilibrium of 
an element Ar cut out of the structure as indicated in 
figure 4, taking advantage of the fact that the shear 
stress is zero at the center line. The shear stress in 
the first panel next to the flange, which is the most 
important one for design purposes, will be obtained 
automatically as part of the solution of the substitute 
structure if i\\e numerical trial -and -error method of 
solution is used, or by using formula (A-2) from ap- 
pendix B in the case of a constant-section panel. 

Panels with variable cross section. In the case of a 
panel with variable cross section, the panel is divided 
into a convenient number of bays as described in 
reference L. For each hay, the cross-sectional area of 
the substitute longitudinal is computed by using for- 
mula (4). In the computation of Kjb by formula (3), 
the average values in the bay are used for A L} b, t, and 
G e /E. The length L is again the total length of the 
panel (not of the bay). The analysis of the sub- 
stitute structure is made by the trial-and-error method 
described in reference L After this step, the procedure 
is identical with the procedure for constant-section 
panels. 

ANALYSIS OF BEAMS WITH FLAT COVERS 

The analysis of beams with flat covers is so closely 
analogous to the analysis of axially loaded panels that 
no detailed example need be given. The substitute 
structure is found exactly as for an axially loaded 
panel. The resulting beam with a single longitudinal 
is analyzed by the trial-and-error method described in 
reference 1, or by formulas if applicable. The total 
force F L at any section can then be distributed over the 
longitudinals as previously described. 
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ANALYSIS OF BEAMS WITH CAMBERED COVERS 

The cambered beam with a single longitudinal. — The 

basic problem of the beam with a single longitudinal and 
a cambered cover was not treated in reference 1. It 
will now be briefly discussed. 

Figure f> shows an element of length d.v cut out of 
the beam. With the help of this diagram, the funda- 
mental equations of equilibrium can be written exactly 
as in the case of the beam with a Hat cover (reference 
l , equations (3a I and (3b) 

<lF F =S w ^-dS c (5a) 

dF L = dS c (5b) 

The equation that expresses the relation between 
shear stress and longitudinal stresses is slightly more 
complicated than in the ease of a flat cover. The 
ordinary bending theory may be taken to give the 
limiting case of no shear deformation. The deforma- 
tions that determine the shear strain must therefore 
be measured from the plane cross section of the engi- 
neering bending theory as a reference base, resulting 
in the equation 

G 

dr= — -gplX**— o>j») — (07.— VLp)]dx (5c) 

where the subscript P denotes stresses obtained with 
the engineering bending theory, which assumes plane 
sections to remain plane. 

These equations can be used to obtain numerical 
solutions by the triaUand-error method, using finite 
differences A in place of the differentials d. Appendix 
B gives the analytical solution for two cases that corre- 
spond to the solutions given for a beam with a (lat 
cover in reference 1 . 
1 In reference l, equation (3b) ts writ ion incorrectly with a minus sign ahead of dfi,. 



z 




THE CAMBERED BEAM WITH MANY LONGITUDINALS 

In the treatment of the cambered heam with many 
longitudinals shown in figure 6 (a), various degrees of 
refinement are possible. The following method, devised 




Figure 5.— Free-body diagrams and notation for single-stringer beam with cambered 

cover. 

to utilize the method developed for axially loaded panels 
and for beams with flat covers, is believed to be ade- 
quate for practical purposes. Attention is called again 
to the basic assumption stated previously, i. e., that 
the camber is moderate. 

The analysis is again divided into two steps: the 
calculation of the flange stresses a F along the span by 
means of the substitute structure, and the subsequent 
distribution of the force F L over the longitudinals at 
any station. 

The area of the substitute st linger is calculated by 
equation (4), using for 6 the developed width of the 
cover sheet. The camber of the substitute beam may, 
for practical purposes, be taken as C s =)4 c (fig. 6 (b)). 
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Figuke G.— Cambered-cover beam. 
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The distribution of the force F L cannot be made 
directly by the method used for flat panels. In the 
flat panel, the longitudinal stress is uniform along the 
chord in the limiting case of infinite shear stiffness; in 
the case of finite shear stiffness, the shear strain is 
defined directly by the longitudinal strains. In the 
cambered cover with infinite shear stiffness, the stress 
varies along the chord according to the straight-line law 
of the ordinary bending theory; in the cambered cover 
with finite shear stiffness, the shear strains are defined 
by the differences between the longitudinal strains and 
the corresponding strains of the ordinary bending 
theory as indicated by equation (5c). 

These differences between the cambered and the Mat 
cover may be interpreted as arising from the fact that 
the cambered cover has bending stiffness of its own 
because it has a "beam depth" equal to its camber. 



.70 
.66 



.54 



























-.30 
























-.34 
















































-.36 

-(!' 
-.42 


















/ 








b * 
































/- 


/ 




















/- 
















-.46 








z. 






















y 











































0 1.0 2.0 3.0 

Yb 

Figure 7— Graph for location of resultant force on cover. 



In the single-stringer beam il was not difficult to take 
care of the effect of this bending stiffness mathematically 
by introducing the terms a FP and a Lr into equation (5c). 
In the multistringer beam it is more convenient to 
introduce a physical equivalent, namely, an auxiliary 
system of longitudinal stresses distributed over the 
cover in such a manner as to make the stress uniform 
in the limiting case of infinite shear stillness. In figure 
6 (e) the broken line shows the stresses given by the 
ordinary bending theory, the full line shows (Ik 1 uni- 
form stress, denoted by a U: and the cross-hatched area 
between the two lines indicates the auxiliary stresses 
necessary to achieve the uniform stress distribution. 
The magnitude of the uniform stress is determined by 
the condition that, when the auxiliary stresses act on 
the flange A F and on the longitudinals A L} they must 
not change the bending moment acting at the section, 
i. c.. they must have zero moment about the assumed 
centroidal line of the lower rover. The auxiliary 
stresses will be denoted by a second subscript A placed 



after the first subscript, which denotes the stringer or 
flange where the stress is measured. 

With the auxiliary stresses assumed active, the 
method of finding the distribution of the stresses along 
the chord is analogous to the method used for Hat 
panels and will he shown in detail for a numerical 
example. From the stresses thus calculated, the 
auxiliary stresses are subtracted to obtain the final 
stresses. 

One step not necessary in the analysis of flat panels 
is required for cambered covers. As indicated in figure 
6 (d), it is necessary to locate the resultant force F L * 
acting on the cover (exclusive of the flange) when the 
actual and the auxiliary stresses are acting. The 
vertical location Ah of this resultant determines tin 4 
effective depth of the beam 

h 4 =h w +M ((>) 

when the combined stresses are acting. The exact cal- 
culation of Ah would require a very tedious integration 
involving the stress distribution and the shape of the 
cover, which has to be repeated several times for each 
cross section with slightly differing values of stress 
distribution. For practical purposes, it will therefore 
be advisable to simplify the problem, although there 
will be a slight loss in accuracy, by finding the lateral 
location y L of the resultant and by assuming that Ah is 
determined by the intersection of the line y=y L and 
the straight line joining F and L. as indicated in figure 
6 (d). Under the assumption of moderate camber, y L 
is given by 

fo **V dA 

Vl= 

fo °* dA 

=~(T5-coth Yb \ csch Yb) 

where )' is (he parameter introduced in reference I, 
equation (21), for the purpose of distributing stresses 
chordwise. The value of y L /b is plotted in figure 7 
against Vb lor ready reference. With the proposed 
simplification, the value of Alt is (hen given by 

Ah = c(l-!/ L /b) (7) 

NUMERIC AL EXAMPLE FOR ANALYSIS OF A 
CAMBERED BEAM 

Figure 8 (a) shows the cross section of the beam 
assumed for the sample analysis. The root section will 
be analyzed for a load P of 250 pounds acting at the tip; 
the length /. of the beam is ins inches. M is assumed 
that the effective width of the sheet has been estimated 
and thai t he value A L 0.85 sq. in. includes the effective 

width. 
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The nexl step is to estimate the effective shear 
modulus. If the presence of camber and of shear 
deformation is neglected, the maximum shear stress in 
the sheet will be given by formula (A 8) of appendix 
B as 

_ PA L _ 250X0.85 _ 0 _ Qn „ , . 
Tmax ~ThA~ T ~0.01l5X2Xl.M~ ' /Sq * 

The buckling shear stress of a long dura] plate 0.0115 
inch thick and 1.80 inches wide is 

r cr «=l,960 lb./sq. in. 

The maximum shear stress being only about twice the 
critical stress, the average shear stress is sufficiently 
close to the critical to neglect diagonal-tension effects 
on shear stiffness and to set G e —G or GJE—0AQ. 
Equation (3) then gives 



KJb 



Vox 



2X0.85X9 _ 



0.535 



.0115X108 2 X0.40 
Inserting this value in equation (4) gives 

A LS =0.85 ^||i = 0.892 sq. in. 

The cross section of the substitute beam is shown in 
figure 8 (b). Since the substitute beam is of uniform 
cross section, it can be solved analytically. Formula 
(A-14) gives 



0.40X0.0115 
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\ 0.80 ^0.892/ 



0.892, 

#=0.0378 KL=4:.08 
Formula (A-12) then gives for z=108 inches 

27, 000X0.94 ^ 1.94X0.892 / 1X0.9991 
? 9.20 L 0.94X0.80 V + 3/T08 J 

o- F =4,830 lb./sq. in. 

This computation completes the first step and the 
substitute beam is discarded. 

The next step is to calculate the stresses in the actual 
beam by the ordinary bending theory: 

^ =27,0 S 0,94gsi2 ' 630 lb - M - in - 

^ = 2 L 0^4 =8220lb/sqin 

Figure 8 (c) shows the chordwise distribution of the 
stresses according to the ordinary bending theory as well 
as the auxiliary stresses. In order to show that the 
auxiliary stresses indicated by figure 8 (c) fulfill the 
requirements, a check on their total moment is made. 

1,955X3 X0.80 = 4,690 
837X3.4X0.189 = 538 

-281X3.8X0.189 = -202 
-1,399X4.2X0.189 =-1,111 
-2,517X4.6X0.189 --2,188 
-3,635X5.0X0.0945 = -1,717 
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The moment is zero with a negligible error. The flange 
stress used for calculating the chordwise stress distribu- 
tion is therefore 

^f* = ^+^=4,830 + 1,955 = 6,785 lb./sq. in. 

The moment furnished by the flange is 

M F *=6,785X0.80X3.00 = 16,280 in.-lb. 




c-2 



A LS -0.892 

',00— "\ 




8,000 



& 6.000 




N (d) 

Figure 8— Cambered-cover beam for sample analysis. 

The moment to be furnished by the cover longitudinals 
is therefore 

Mz,* = 27,000- 16,280- 10,720 in.-lb. 
Assuming /< c =3.77 inches, the force F L * becomes 

/^* = 1 -M20 =2 ,840 lb. 

The average stress is therefore 

* 2,840 Q QA , 

and the Patio 



GL av 3,345 n 

"^-6785 

With this value as abscissa, read from figure 18 of 
reference 1 

T7> = 1.94 



80^93— 38- 
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and, with this value, 

7^/6 = 0.613 

from figure 7 so that 

AA=2(1— 0.613) =0.77 in. 
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The stress at the center line is given by 



O 2 4 6 8 10 

oc/L 

Figure 9.— Stresses in axially loaded panel (experimental data from reference 2). 

and 

fc e =3.00+0.77==3.77 in. 

which agrees with the assumed value. If it did no1 
agree, a new trial would have to be made. 



<?CL 



(i,7sr, 



= 1,910 



cosh Yb 3.551" 

From this stress the actual stress is obtained by sub- 
tracting the auxiliary stress 

<r ei *= 1,910— (— 3,635) =5,545 lb./sq. in. 

Table I shows the calculation of the stresses a-* in the 
other stringers by the formula 

%*=oW |t cosh Yy 
and of the final stresses; figure 8 (d) shows graphically 
the final stress distribution. 

TABLE I 



Stringer 


m 


Yy 


cosh Yy 


a* 

(lb./sq. in.) 


OA 

(lb./sq. in.) 


(lb./sq. in.) 


Center line 


0.0 


0 


1.000 


1.910 


-3, 635 


5,545 


1 


.2 


.388 


1. 076 


2,055 


-2, 517 


4, 572 


2 


.4 


.776 


1. 316 


2, 514 


-1,399 


3,913 


3 


.6 


1. 164 


1.758 


3. 360 


-281 


3,641 


4 


.8 


1.552 


2. ir.ii 


4,710 


837 


3, 873 


Flange. 


1.0 


1.94 


3. 551 


6, 785 


1, 955 


4, 830 



EXPERIMENTAL STUDIES 

AXIALLY LOADED PANELS 

Experimental results for a panel loaded in compres- 
sion are described in reference 2. This panel was shown 
in figure 3 and served as a numerical example for the 
proposed method of analysis. The results of the 
analysis as well as the experimental results are shown 
in figure 9. 

GENERAL REMARKS ON ANALYSIS OF BEAM TESTS 

In the analysis of beam tests, some difficulty is met 
in establishing the idealized section. It is easy to 
define locations for the longitudinals but fixing the 
location, and particularly the size, of the Manges pre- 
sents diMiculties, because part of the shear web must 
be considered as furnishing a ^contribution to the 
idealized flange. 

In order to reduce arbitrariness to a minimum, the 
following procedure was adopted for all beam analyses. 
First, the centroidal axis and the geometric moment of 
inertia of the cross section in question were computed. 
If the sheet was considered to be only partly effective 
in carrying normal stresses, the proper effective width 
was used in these computations. Next, the locations of 
the idealized Manges were fixed. On that side of the 
beam where the shear deformation was being calculated, 
the flange was assumed to be in the plane of the cover 
sheet. On the other side, which was without cover 
except in one case, the flange was assumed to be located 
at its estimated centroid. The cross-sectional areas of 
the tw r o idealized flanges (tension and compression) were 
then computed from the conditions that the idealized 
section must have the same centroidal axis and the same 
moment of inertia as the actual section. For this 
idealized section, the analysis was then made by using 
the previously described methods. 
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GENERAL REMARKS ON N. A. C. A. BEAM TESTS 

The available published test data were not sufficient 
for ;m adequate check of the theory developed. A 
number of beams were therefore tested by the N. A. C. 
A. The first of these beams was discussed in reference 



decreased again in steps of 100 pounds. The slope of 
the straight line through the test points was used to 
determine the stress at P=250 pounds, which will be 
shown in the later figures. 



2xj 1 7 ST, 



3' 2.02 lb 
17 ST beam 



,0.0115 17ST 



,T\no 2xj 17ST 



4.5 



Rivet 
spacing 



Rivet 
Spocinq 



^2xj 17 ST 



r- 4.5 ~ 4.5 > 
A'O.85 


A F -0.85 


r 9 


ij AL-Q85. 6 


1 


^-0.0115 


Ac 


A 


Ge/E-0.288 


(b) 


(c) 





1907 



S 



1 I i 1 1 1 1 

x Experimental point 's I 

8 Coinciding experimental points 

— Co/cu/ated, first approximation 

— * , second 



Mc/I 

-Younger's solution 



(a) Actual section, 
(b) Idealized section. (c) Substitute section. 

Figure 10— Cross section of N. A. C. A. beam 2. 



1 ; the following ones, designated as N. A. C. A. beams 
2, 3, and 4, will be discussed in this paper. 

In all N. A. C. A. beams, measurements were made on 
the tension side of the beam in order to eliminate 




Figure 11.— N. A. C. A. beam 2 under test. 

erroneous strain readings caused by local buckling of the 
stringers. Furthermore, flat strips could be used for 
stringers, making it possible to take strain readings very 
close to the sheet. 

The load 2P was increased from 0 to 500 pounds (in 
the first series of tests) in steps of 50 pounds and 




10,000 



8,000 



6,000 



80 60 40 

Distance from tip. in 

Figure 12— Stresses in N. A. C. A. beam 2 for P=250 lb. #=10.4X10°. 
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Readings were taken across the entire section of the 
beam and on both sides of the stringers. Each point 
representing a flange stress in the figures is therefore the 
average of two slopes, and each point representing a 
stringer stress is the average of four slopes excepting 
figures that show the stress distribution along the 
entire chord. 

A slight departure was made from the described 
method of analysis in the case of N. A. C. A. beams 2 
and 3. The uniform distribution of A L along the chord 
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Figure 13 —Cross section of Galcit beams, a, upper cap angle. Area= 0.354 sq. in. 
b, beam web. *=0.051 in. c. lower cap angle. Area=0.221 sq. in. d, stiiTener 
angle. Area=0.044 sq. in. e, attaching strip. Area=0.082 sq. in. f, cover sheet. 
fi=0.025 in.; f 2=0.050 in. 

is not very well approximated in these beams, A L con- 
sisting of only two stringers. Consequently, equation 
(4) was not used. The two stringers constituting A L 
were treated individually on the basis of equations (1) 
to (3). This departure also accounts for the fact that, 
for beam 3, the substitute camber is not taken as one- 
half the actual camber, as recommended for practical 
cases with many stringers. For comparison with the 



experimental flange stresses, the stresses calculated for 
the idealized flanges of the N. A. C. A. beams were 
corrected to the outside fiber stresses on the assumption 
that plane sections remain plane. For the purpose of 
calculating the shear deformation, the width of the 
sheet was taken between rivet rows for N. A. C. A. 
beams 2 and 3. 

TESTS ON BEAMS WITH FLAT COVERS 

N. A. C. A. beam 2. — N. A. C. A. beam 2 was similar 
in design to beam 1 described in reference 1. The 
cross sections of the beam are shown in figure 10. The 
bulkheads, not shown in this figure, were similar to 
those on beam 1 and were spaced to make the bays 
about square. The length L of the beam was 108 
inches. Figure 11 shows the beam under test and 
figure 12 shows the results of the tests and of the 
calculations. 

Galcit test beams. — Figure 13 shows the cross section 
of a type of beam tested at the California Institute of 
Technology (reference 3) under a pure bending moment. 
Figure 14 (a) shows the experimental and calculated 
results for the beam with t— 0.025 inch and figure 14 (b) 
shows the results for the beam with t= 0.050 inch. 

Schnadel's ship model. — Figure 15 shows the cross 
section and the side view of a ship model tested by 
Schnadel (reference 4). The model was built of steel. 
Measurements were taken only on the outside of the 
compression cover (corresponding to the deck of the 
vessel) over one quadrant of the beam. Figure 1G 
shows the experimental results and the results calculated 
by the method presented in this paper. 




eo o ?o 

Distance from root, in. 
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(a) Beam 1; *=0.025 in. (b) Beam 2; <=0.050 in. 

Figure 14.— Stress distribution in Galcit beams. Experimental data from reference 3; A/b=120,000 in.-lb.; G,/E assumed 0.25 (diagonal tension). 
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TESTS ON BEAMS WITH CAMBERED COVER 

N. A. C. A. beam 3. — Figure 17 shows the cross sec- 
tions of a cambered beam obtained by inserting cam- 
bered bulkheads into N. A. C. A. beam 2. Figure 18 




: 40 =+ 



in 



w (b) 

(a) Cross section. 



80- 



T 



,5/ of ions 



® 

(b) Loading diagram. 



Figure 15.— Cross section and general diagram of Schnadel's ship model. Dimen- 
sions are in cm and loads in kg. 

-hows the beam under test. Figure 19 is a view of the 
inside of the beam, showing intermediate bulkheads 
that were added for tests at high loads to reduce sagging 
of the stringers between the main bulkheads. This 
sagging is proportional to the square of the stresses 
and consequently may become important at high 



stations for three different loads. Two facts are evident 
from an inspection of this figure: The differences be- 
tween the actual stresses and the stresses of the ordinary 
bending theory increase as the root is approached and 




3" 2. 02 lb. 
17ST beam 



N20//5 \7ST 
Vivo 2*i 17ST- 

acmq < 



Rivet 
spacing 



-2yj- \7ST 



fa) 




G*/E- 0.888 
W (c) 

(a) Actual section, 
(b) Idealized section. (c) Substitute section. 

Figure 17— Cross section o! N. A. C. A. beam 3. 

also as the load increases, because the shearing stiffness 
decreases with increase in load. 

N. A. C. A. beam 4. N. A. C. A. beam 4 was tapered 
in plan form, in depth, and in stringer area as shown in 
figure 22. Figure 23 shows the calculated and experi- 
mental stresses in the flange. The experimental 
stresses shown in this figure are based on measurements 
taken on the outside of the flange but are corrected to 
the top edge of the web. 
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Figure 16 — Stress distribution in Schnadel's ship model (experimental data from reference 4). 



stresses, but it requires attention only in the case of 
shallow beams. Figure 20 shows experimental and 
calculated stresses in this beam at P=250 pounds. 

A shorter series of measurements was made on beam 3 
at higher loads. Figure 21 shows the stresses at four 



Figure 24 shows the chordwise stress distribution at 
the station z=91.4 inches. The experimental stresses 
shown are not the stresses measured on each stringer 
but are weighted averages of the stresses measured on 
each stringer and on the skin adjacent to the stringer on 
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each side. It was found fchat the skin stresses were 
consistently higher than the stringer stresses; near the 
root the difference was as much as 20 percent, but the 
difference decreased (roughly proportionally) with 
distance from the root. Since the ratio of stringer area 
to sheet area was more than 4:1, the weighted average 



ences between stringer stresses and skin stresses, pre- 
viously mentioned, are considered. 
For the Galcit beams, the agreement is somewhat 




Figure 18.— N. A. C. A. beam 3— view of closed side. 

stress never differed by more than 5 percent from the 
stringer stress proper. 

DISCUSSION OF RESULTS 

COMPARISON BETWEEN PROPOSED METHOD OF ANALYSIS AND 
EXPERIMENTAL RESULTS 

The agreement between experiment and calculation is 
good for the axially loaded panel (fig. 9). For N. A. 




Figure 19.— N. A. C. A. beam 3— view of open side. 

C. A. beams 2 and 3, the agreement is good except for 
the root region of the center stringer in beam 3 (figs. 12 
and 20). 

For X. A. C. A. beam 4, the agreement is reasonably 
good for the flange stresses (fig. 23). For the stringer 
stresses, which are shown only for the root station in 
figure 24, the agreement may be considered fair, if the 
differences between the two test series and the differ- 



I r — i — 

x E xperimental points 
H Coinciding experimental points 
Co/cutoted, with sheor deformation 
by Mc/J 




8,000 



60 40 
Distance from tip f in. 

Figure 20— Stresses in N. A. C. A. beam 3 for P=250 lb. #=10.4X10«. 

poor for the beam with the thin cover (fig. 14 (a)) but is 
quite good for the beam with the thick cover (fig. 14 
(b)>. 

For SchnadePs ship model, the agreement is fair at 
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some stations and very poor at others. Study of the 
test report shows that the accuracy of the test was, for 
a number of reasons, far below the accuracy of all 
other tests analyzed in the present paper. This con- 
clusion is borne out by inspection of the results in 
figure 16. Note, for instance, at station 2 and particu- 
larly at station 3, that all experimental stresses are 



22,000 



20,000 



/ 6,000 



16,000 



/4,00Q 



t>/ '2,000 



beams under load show s thai the spanwise variation in 
the condition of the shed is indeed small; it should be 
borne in mind thai relatively large variations in shear 
stillness influence the stringer stresses but little, as 
shown in reference 1. The chordwise variation, how- 
ever, is marked, the outer panels being buckled while 
the inner panels are not. This variation was taken 
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(a) .r= 103.5 in. (b) x»97.5 In. (c) j=f>7.5 in. (<)) 1=59.5 in. 

Figure 21.— Chord wise stress distribution in N. A. C. A. beam 3 at four stations for P=260, 500, and 760 lb. 
^fc/I shown for P=25() and 750 pounds. 



considerably higher than the calculated ones; hence 
the summation of the interna] moments would be much 
larger than the external moment. The test was in- 
cluded in the analysis because it is the only available 
complete test on the limiting case where stringers and 
sheet are merged into a single unit, a plate. 

In all beam analyses made for the present paper, 
over-all average values of effective shear stiffness were 
used. A glance at the photographs of the N. A. C. A. 



into account approximately by using a weighted average 
value of Og, and this procedure may be responsible for 
some of the discrepancies between test and calculation. 
Theoretically, it might be possible to take this variation 
into account more exactly, but there appears to be 
little justification to do so when the proposed simplified 
method of analysis is used. In practical design, large 
chordwise variations of shear stiffness should be 
avoided by using heavier skin near the flanges. 
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If the far-reaching simplifications involved in the 
theory are considered as well as the difficulties of strain- 
gage testing of sheet-metal structures, the agreement 
between experiments and analysis is, on the whole, 
fairly satisfactory. Although the analysis does not 
give a perfect picture of details, it does appear to give 
a substantially correct picture for the stresses most 
important in design work. 

To persons unacquainted with strain-gage testing, 
the discrepancies between tests and calculations might 
appear to be rather large. It should be pointed out, 
how r ever, that strain-gage tests of conventional types 
of structures, such as trusses and plate girders, fre- 
quently show discrepancies fully as large or larger. 
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Figure 22— Cross sections of N. A. C. A. beam 4. 

COMPARISON BETWEEN PROPOSED METHOD OF ANALYSIS AND 
YOUNGER'S SOLUTION 

The proposed method of analysis is based on the 
-nine simplified physical concepts as Younger's method 
deference J). Younger's solution is mathematically 
more rigorous, but it applies only to a beam of constant 
section with a cosine-wave bending moment. For 
practical shapes of bonding-moment curves, it is neces- 
sary to superpose a number of cosine terms. 

Comparisons for the case of a concentrated load ap- 
plied at the tip show that the substitute-structure 
method of analysis gives flange stresses at the root that 
are as much as 15 percent higher than the stresses cal- 
culated by superposing four cosine terms. Judging by 
the magnitude of successive terms, four were considered 
a sufficient number to give the desired accuracy. The 



stresses in the longitudinals given by the substitute- 
structure method are correspondingly lower than those 
given by Younger's formula. Comparison with experi- 
ments for two cases (figs. 9 and 12) shows that Younger's 
extended formula is in very much poorer agreement 
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Figure 23.— Flange stresses in N. A. C. A. beam 4 at P=250 lb. JE=10.4X10 6 . 

with the experiments than the substitute-structure 
method. This fact is somewhat surprising, and the 
question arises as to what might be the possible reasons 
for the poor agreement. 

If a diagonal-tension field forms on the sheet, the 
shear between flange and longitudinals will not be 
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Figure 24.— Stresses at first station in N. A. C. A. beam 4. 

transmitted at right angles to the axis oi the beam but, 
theoretically, at 45° angles. The theory may therefore 
be expected to give reasonably accurate results only if 
the bending moment does not change too much over a 
spanwise distance equal to the width of the beam. 
Obviously, this condition is not fulfilled by the higher 
cosine terms after the first one, so that their physical 
significance may be seriously questioned. 
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The evidence presented by the axially loaded panel 
(fig. 9), which did not buckle appreciably at P= 1,000 
pounds, appears to indicate that even for a shear-resist- 
ant panel the superposition of cosine terms does not al- 
ways yield a sufficiently close approximation to the phys- 
ical facts. If this defect always exists, then any method 
based on the same fundamental physical concepts and 
relying on trigonometric series will be unreliable. 

It might be mentioned in passing that the theoreti- 
cal treatment given by Schnadel in reference 4 and in 
several other papers is of little interest for aeronautical 
structures, because it applies only to an isotropic plate 
where the shear stiffness is fixed by the theoretical 
relation 

^"~2(1 + M ) 

The results therefore contain no provision to take into 
account reduced values of G e or plates stiffened by 
st lingers. 

THE INFLUENCE OF RIBS 

Ribs or bulkheads influence the stresses in the beam 
cover in two ways. By virtue of their axial stiffness, 
they help to carry transverse stresses in the cover. 
This function is unimportant if the sheet does not 
buckle into a diagonal-tension field, but it is, of course, 
of paramount importance if a diagonal-tension field 
forms. Because the rib flanges have bending stiffness 
in the plane of the cover, they also tend to reduce the 
shear deformation. It was pointed out in reference 1 
that this effect can be calculated for a single-stringer 
beam and it was stated that in practical cases the effect 
is very small. This conclusion, drawn from the calcula- 
tions, has been confirmed by tests of N. A. C. A. beam 2. 
It should be noted, of course, that these remarks apply 
only if the basic requirement of very moderate camber 
is fulfilled. 

N. A. C. A. beam 2 was tested first with all longitu- 
dinals sliding freely over the ribs and held against the 
ribs only by their ow n tension. A second strain survey 
was then made of the beam after connecting the longi- 
tudinals with the ribs by taper pins. The ribs were 
very heavy steel channels, as shown on the drawings 
and photographs, but their only effect was to smooth 
out a few minor irregularities in the stress-distribution 
plots. An extremely heavy tip rib was then added; 
this rib reduced the stress in the flange about 6 percent. 
Calculation indicated, however, that an equivalent 
amount of material used to thicken the skin would have 
resulted in increasing the skin thickness by about 500 
percent over the entire span and would have reduced 
the stress in the flange by about 33 percent. 

A brief inspection of figure 25 is sufficient to show why 
the rib is quite ineffective. Figure 25 (a I shows the tip 
rib acted upon by the longitudinal. In figure 25 (b) it 
was assumed that the material contained in the tip rib 
is spread out some distance along the span. It is obvi- 
ous that this change results in a much stiffer cross beam. 



All tests of N. A. C. A. beam 3 were made without 
connections between longitudinals and ribs. On N. A. 
C. A. beam 4, which had bulkheads of normal size, the 
longitudinals were riveted to the bulkhead flange. 

THE EFFECTIVE SHEAR MODULUS 

The effective shear modulus of a thin sheet framed 
by rigid edge members is equal to the shear modulus of 
the materia] as long as the shear stress is lower than the 
critical or buckling st ress. If the stress is increased be- 
yond this value, diagonal-tension folds begin to form 
and grow. The effective shear modulus gradually de- 
creases, approaching asymptotically the value G e ~ 
b /$ G. The nature of this transition was investigated 
experiment ally by Lahde and Wagner (reference ">). 

In practical structures, the edge members are not 
rigid; they have a finite axial stiffness and a finite 
bending stiffness. The influence of these stiffnesses 
has been treated analytically by Wagner in his original 
theory for the ease of a fully developed diagonal-tension 
field. The influence of edge members with finite stiff- 
ness on the characteristics of a thin sheet in the transi- 




(a) 



^^^^^^^^^ 

(b) 



Figure 2.5. 



tion zone between shear-resistant sheet and diagonal- 
tension field has not been investigated to date. At kin 
offers a method of estimating the characteristics of a 
diagonal-tension beam by making tests on square panels 
(reference 6). Although this idea is fundamentally 
sound, At kin's analysis is open to a serious objection. 
He claims that the deflection 5 of a test panel can always 
be represented as a straight-line function of the load 

/'. and he sets 

8=kP 

which means, in effect, that Atkin's method takes into 
account only the finite stiffness of the edge members. 
It disregards the gradual transition from G e =G to 
G e = %G in a rigidly framed sheet. In many practical 
cases, where the critical stress is not exceeded more than 
three or four times, the second factor is probably far 
more important than the first. 

The tests described in references 2 and 3 were evalu- 
ated by their authors to give values of effective shear 
stiffness. These analyses have been questioned in a 
later paper (reference 7), chiefly because the values 
obtained were much lower than the theoretical values for 
the pure diagonal-tension Held with rigid edge members. 
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A critical examination shows that in all those 4 analyses 
the shear stillness has heen obtained by taking the 
differences of slopes of two experimental curves. This 
method is extremely sensitive to slight experimental 
errors. Unfortunately, experimental errors in strain- 
gage tests of sheet-metal structures are quite large, and 
the stresses are, furthermore, quite insensitive to 
changes in shear stiffness. The results obtained by such 
a slope method are therefore very questionable, and in 
some cases it is possible to change the calculated values 
of the shearing stiffness several hundred percent by 
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varying, for instance, the effective width within its 
possible limits. 

In view of these circumstances, it appears more 
advisable to analyze 4 the tests in such a manner that 
over-all average values for the shear stiffness are ob- 
tained by utilizing the ordinates of experimental curves 
instead of the slopes of these curves. The procedure 
would be to calculate the stresses under several assump- 
tions for the shear stiffness and to find the stress curve 
thai gives the best agreement with the test results. 
Several examples of such a procedure are given in refer- 
ence 1. 

The high-load tests of N. A. C. A. beam 3 were 
analyzed in a similar manner, I nfortunately, the 
limited number of strain gages necessitated repeat load- 



ings; during these repeat tests, changes occurred in 
parts of the beam that prevented a definite analysis. It 
was estimated that, at P=900 pounds, the effective 
shear modulus was G e =0.5G, but no definite estimate 
could be made for higher loads. The trouble may have 
heen partly that the flange was no longer obeying 
Hookc's law at the gage station, the stress being over 
30,000 pounds per square inch. 

It should be noted that the effective modulus was 
well below the theoretical value G e =%G for rigid edge 
members, in spite of the fact that the edge members 
were much stiffer than they would be in actual construc- 
tion and that tension was superimposed on the shear in 
the skin. 

The opinion is occasionally heard that the shear 
modulus of corrugated sheet is appreciably less than 
that of Hat sheet. There appears to be no published 
information to support such an opinion. The analysis 
of torque tests of box beams with corrugated covers 
(reference 8) leads to the conclusion that up to shear 
stresses of around 3,000 pounds per square inch the 
shear modulus of corrugated sheet is equal to the 
modulus of the material. Small deviations of 5 to 10 
percent, which occur in such tests, can probably be 
attributed to inefficiency of the joints in the built-up 
boxes because they have 4 been found in practically all 
torque tests. Kbner, who has an exceptionally broad 
background of experience in tested stressed-skin struc- 
tures, states iti reference 9 that the sheer stiffness of 
corrugated sheet remains unchanged up to the point of 
failure. It is necessary, of course, to make proper 
allowance for the difference between developed width 
and projected width of a corrugated panel when com- 
puting shear deformations. 

APPLICATION OF THEORY TO FUSELAGES 

The theory in this paper was developed for the ex- 
press purpose of furnishing means for analyzing wing 
beams or other beams with very moderate camber. It 
is of interest, of course, -to gain some idea of how well 
the theory applies to beams with large camber, such as 
fuselages. A fuselage test that came to the attention of 
the author after the investigation was finished will 
therefore be included. 

The details of the test may he found in reference L0. 
The most important data are given in figure 2(5. The 
shell represents a fuselage with symmetrical cut-outs, 
and the bending moment is introduced in the form of 
concentrated forces at the longerons. Between frame a 
and the end, the shell was fixed to a test jig by a heavy 
steel ring. 

The part of the shell between the longerons and the 
neutral axis was considered as "shear web" and the 
remainder as "cover." The analysis was made by 
formula (A-16). Local corrections to the computed 
stresses were made between frames e and g, because the 
sheet thickness was 0.08 centimeter between e and 
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f, and 0.10 centimeter between f and g resulting, to- 
gether with the cut-out, in some changes in effective 
areas in this region. 

The only variation from the standard procedure out- 
lined in this paper was the use of a somewhat more 
rational method of determining c& than simply assum- 
ing c s =%c. The substitute camber was determined 
by the condition. 




8 f e d c b a 

Frames 

Figure 27.— Experimental and calculated stress distribution in fuselage with sym- 
metrical cut-outs (experimental data from reference 10). 



That is, if the longitudinals ;ire concentrated at the z 
location defined by c s , the moment of inertia must be 
the same as in the actual section. 

Figure 27 shows the experimental and the calculated 
results. 

CONCLUDING REMARKS 

Large shear deformations are probably always ac- 
companied by loss of structural efficiency; efficient 
design therefore calls for utilization of all available 
means for reducing the shear deformation. In a com- 
bination consisting of skin stiffened by individual 
stringers, the stringers furnish no contribution to the 
shear stiffness. In a combination of flat skin with 
corrugated skin, however, all the material carries longi- 
tudinal stresses as well as shearing stresses; such a 



combination probably represents, therefore, a close 
approach to the best possible efficiency from consider- 
ations of uniform stress distribution. It must be 
remembered, too, that the shear stiffness of flat sheet 
is very adversely affected if it is thin enough to buckle 
into diagonal-tension folds, a condition that does not 
develop in corrugated sheet. 

For sheet with individual stringers, experimental 
studies on individual panels have usually led to the 
conclusion that the best efficiency is obtained by mak- 
ing the skin as thin as possible, consistent with practical 
considerations. If the shear deformation in the actual 
structure is taken into account, it becomes evident that 
this conclusion will often require serious modification. 
It might be worth while in some cases to investigate the 
effect of thickening the skin near the wing tip, where the 
shear deformations are largest and therefore easiest to 
decrease. It might be pointed out that, once an ade- 
quate tip rib is provided, shear deformation can be 
reduced more efficiently by increasing the skin thickness, 
especially near the tip, than by attempting to increase 
the (horizontal) bending stiffness of the tip rib. 

A final word of warning should be given. A method 
of stress analysis such as the method described in this 
paper deals only with the stress distribution before 
failure occurs. If the maximum stress for a given load 
is varied by changing the design of the structure, then 
the failing stresses may change, too, so that the maxi- 
mum stress is not the sole criterion for the efficiency of 
the structure. For example, if the skin is made very 
heavy with relation to the stringers, then buckling of 
the skin may induce premature failure of the stringers. 
Thus far, no mathematical analysis of this problem has 
been published; test results must be used. The subject 
of allowable stresses is beyond the scope of this paper, 
but is mentioned herein as a warning against drawing 
hasty conclusions. 
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APPENDIX A 



LIST OF SYMBOLS 

Ay cross-sectional area (sq. in.). 
Ej Young's modulus (lb./sq. in.). 

F, internal force (lb.). 

G, shear modulus (lb./sq. in.). 
/, geometric moment of inertia. 
K, constant. 

L, length of panel or beam (in.). 

.1/, bending moment (in. -lb.). 

/\ externa] load (lb.). 

S, shear force (lb.). 

bj half-width of beam or panel (in.). 

C, camber of cover (in.). 

//, depth of beam (in.). 

t, thickness of cover sheet (in.). 

running load (lb. /in.), 
a;, dist ance along center line. 
y, distance from center line. 
.:, distance from centroidal axis of cross section. 
18 



a, direct (normal) stress (lb./sq. in.), 
r, shear stress (lb./sq. in.). 

*, denotes condition where actual and auxiliary 
stresses are superposed. 

Subscripts have the following significance: 

A, auxiliary. 

C, cover sheet. 

CLy center line. 

F, flange. 

L, longitudinal. 

P, theoretical values assuming that plane sections 
remain plane. 
Sy substitute. 
T, total. 
TJy uniform. 
Wy shear web. 
a, applied. 
e, effective. 
0, root section. 



APPENDIX B 



ANALYTICAL SOLUTIONS FOR STRUCTURES WITH A 
SINGLE LONGITUDINAL 

GENERAL REMARKS 

The sign conventions of reference 1 are retained. 
Stresses in stringers are positive when tensile. Shear 
stresses in the cover sheet are positive when caused by 
positive stresses in the flange. F. Shear stresses in the 
shear web are positive when causing positive stresses in 
the flange F. 

The figures show, first, the half structures and, sec- 
ond, the two possible cases of making symmetrical 
structures out of these half structures. The formulas 
should be applied only to such symmetrical structures. 
Theoretically, the formulas also apply to the half struc- 
tures if the forces T are applied at the stiff transverse 
member a1 the tip. This procedure would involve the 
assumption that the stringers were infinitely stiff in 
bending; it is therefore believed that the application of 
the formulas to the half structures might easily lead to 
very serious errors. 

Some of the formulas have already been given in ref- 
erence 1. They are repeated here for convenience and 
are written in a slightly different form to bring out more 
clearly the correction factor that must be applied to the 



(a) 



fb) 
Figure 28. 
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ordinary bending theory in order to take shear defor- 
mation into account. When the shear stiffness ap- 
proaches infinity, this correction factor approaches zero. 



I— THE AXI ALLY LOADED PANEL 



(a) The longitudinal built in at the root (fig. 28).— 
For the case of an axially loaded panel with the longi- 



tudinal built in at the roof, the following formulas are 

obtained: 

Let 
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(a) 



and 

Then 
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(b) 
Figure 29. 



A T = A F +A L 
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T ~A F EbK cosh KL 



'///////////A/////////////* 
U 



(c) 



A L cosh Kx 
A F cosh KL 



_ P( cosh IS A 
aL ~A T \ cosh KL J 
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(A-2) 
(A-3) 
(A-4) 



(b) The longitudinal not built in at the root (fig. 
29). — The easiest way to treat the case of the longi- 
tudinal not built in at the root is to take advantage of 
the symmetry of the structure. When the origin is 
taken at the middle of the length L, this case is reduced 
to case I (a). 

II— THE BEAM WITH FLAT COVER 

The formulas for the beam with flat cover apply to 
two cases: beams in which the depth h is constant along 
the span, if a concentrated load P is applied at the tip; 
and beams in which the depth // tapers linearly to zero 
at the tip, if the loading w per (oot run is uniform along 
the span. In the case of uniform loading, wL/2h Q is 
substituted for Pjh in the formulas for shear stress. 
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4. PROPELLER SYMBOLS 



D, Diameter 

p, Geometric pitch 

p/D, Pitch ratio 

V, Inflow velocity 

V 8t Slipstream velocity 

T 

T, Thrust, absolute coefficient C T — — rru 

Q 9 Torque, absolute coefficient C Q = p ^jys 



p, 

c„ 

v, 
n, 



Power, absolute coefficient C P = 



4 



Pri 2 



Speed-power coefficient = 

Efficiency 
Revolutions per second, r.p.s. 

Effective helix angle =t an" 



5. NUMERICAL RELATIONS 



1 hp. = 76.04 kg-m/s=550 ft-lb./sec. 
1 metric horsepower= 1.0132 hp. 
1 m.p.h. =0.4470 m.p.s. 
1 m.p.s. =2.2369 m.p.h. 



1 lb. = 0.4536 kg. 

1 kg=2.2046 lb. 

1 mi. = l,609.35 m=5,280 ft. 

1 m=3.2808 ft. 



